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Abstract 

The spin structure function of the proton is studied in a two component 
framework, where the perturbative evolution of parton distributions and nonper- 
turbative vector meson dominance model are used. We predict the asymmetric 
behavior at small x from lower to higher Q^. We find that the contribution of 
the large gluon helicity dominates at x > 10 “^ but mixed with nonperturba- 
tive component which complicates the asymptomatic behavior of aX x < 10 “^. 
The results are compatible with the data including the HERA early estimations 
and COMPASS new results. The predicted strong Q^- and x-dependence of < 7 ^ at 
0.01 < < 3GeV^ and x < 0.1 can be checked on the next Electron-Ion Collider. 
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1 Introduction 


Recently, COMPASS experiment at CERN collected a large number of events of polar¬ 
ized inelastic scattering off the protons with very small values of Bjorken scaling variable 
X [1]. The preliminary analysis of these data combining with the previous experiments 
[2], showed non zero and positive asymmetries of the structure function gf. In these hxed 
target experiments the low values of x are almost reached by lowering the values of Q^. 
The knowledge of the nucleon spin structure function gi{x,Q^) at low and small x 
is particulary interesting, since it is not only an important information to resolve the 
’’proton spin crisis”, but also provides us with a good place to study the transition from 
the perturbative research to the nonperturbative description of the proton structure. 

In this work we try to study the behavior of gf at small x but in the full range. As 
we know that the structure functions of the nucleon are mainly constructed of the parton 
distributions at > lGeV‘^, while the non-perturbative contributions to the structure 
functions become unneglectable at 1 GeV^. A key question is what components 

construct the spin structure functions of the proton at such low Particularly, do the 
parton distributions and their pQCD evolution still play a role or not? For answering these 
questions, we introduce the application of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) equation [3] with the parton recombination corrections at low in detail. 
These corrections have been derived both for the polarized and unpolarized parton distri¬ 
butions in our previous works [4,5]. We point out that the isolation of the contributions of 
the vector meson is necessary for keeping the factorization schema of the polarized parton 
distributions at low Q^. We hnd two different asymptotic behaviors of at a; < 10“^: 
nonperturbaive behavior ~ x~^ at < IGeV"^ and perturbative drop at > ?>GeV^. 
We predict the translation of gi at small x from lower to higher . The results 
are compatible with the data including the early HERA estimations and GOMPASS new 
results. We point out that the measurements at different x with different values of 
in the hxed target experiments mix the complicated asymptomatic behavior of g^. The 
predicted strong Q'^- and x-dependence of gi at 0.01 < < 3GeV^ and a; < 0.1 due to 

the mixture of nonperturbative vector meson interactions and the QGD evolution of the 
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parton distributions can be checked on the next Electron-Ion Collider (EIC). 

The organization of this paper is as follows. In Sec. 2 we discuss the applications of the 
polarized parton distributions at low based on the generalized leading order approx¬ 
imation. In Sec.3 we summarize the contributions of the polarized parton distributions 
to the spin structure function gf of the proton, which have been hxed by our previous 
work. The contributions of the vector meson to Qi are discussed in Sec. 4. We present 
our predictions of Qi and the comparisons with the data in Sec. 5. The discussions and 
summary are given in Sec. 6. 
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2 A general consideration of the nucleon structure 

function at low 

In the researches of the nucleon structure functions at the full kinematic region, an 
argued question is whether the parton distributions and their perturbative QCD evolution 
can (even partly) be applied to the low range or the parton concept is suddenly invalid 
at a critical value of < IGeV^^ 

Let us begin from the parton model for the spin-dependent distribution, which is 
written based on the Collins-Soper-Sterman (CSS) factorization schema [6] at the collinear 
approximation and the twist-2 level, 

9i{x,Q^)= f —J2Cqix/y,Q^/nF)Sqiy,^F), ( 2 . 1 ) 

Jo y q 

which breaks up the spin structure function into two factors associated with perturbative 
short-distance functions Ca and nonperturbative polarized parton distributions 5q at the 
factorization scale y,F- 

Taking the lowest order of Cq 

C,(x/y, - 1)^(0 - yp) + 0(a.) + 0(l/<3), (2.2) 

0{as) and 0(1/ Q) are the QCD radiative corrections and higher twist contributions. 
Inserting it to Eq. (2.1), we obtain the relation between the spin structure functions and 
the polarized quark distributions 

gi{x, Q^) = X X! ^q[^9{x-, Q^) + Sq{x, Q^)] + 0(as) + 0(1/Q). (2.3) 

^ <7 

According to the renormalization group theory, 

dgpO^ = 0 , ( 2 , 4 ) 

dlnyF 

it gives the DGLAP equation 
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(2.5) 


Q‘^^Sq{x,Q^) = j —J2^Pqq'{x/y,^s{Q^))Sq'{y,Q‘^), 

dQ^ Jo y ^ 

APqq/ denotes the splitting functions. If we consider only the leading order {LO) approx¬ 
imation, we have 


gDGLAP(^^^ Q‘^) = ^J2el [Sq{x, Q 2 ) + Sq{x, Q^)] , (2.6) 

These results are available at > a few Gel^^. 

Now let us consider what will happen to the above results at low In principle, the 
initial parton distributions dq^x^jj?) in Eq. (2.1) are defined at = fj? {Aqcd < h << 
IGeE). Therefore, the parton distributions are non-perturbative essentially. However, the 
factorization in Eq.(2.1) should be modihed at low because of the following reasons: 

(i) The handbag diagram Fig. la is a typical time ordered diagram describing Eq. 
( 2 . 1 ), where the quark propagators connecting with the probe and the target have only 
the forward component, these propagators can be broken as shown in Eq. (2.1) since they 
are on-mass-shell. The corresponding backward quark propagators construct the cat’s 
ear diagram Fig.lb, which are neglected since these backward propagators are absorbed 
by the target in the collinear approximation [7]. However, the contributions of Fig. Ic 
to Eq. (2.1) can not been neglected at low due to the corrections of quark-antiquark 
pair, which interacts with the target as a virtual vector meson if the transverse momentum 
k± ^ Q oi quark pair is not large and conhnement effects are essential. The interference of 
the forward and backward quark propagators in Fig. la and Ic will put these propagators 
to off-mass-shell and breaks the factorization schema. To avoid this event, we use a 
phenomenological vector meson dominance (VMD) model [ 8 ] to ’’isolate” the contributions 
from Fig. Ic. Traditionally, such VMD hypothesis is used to explain the structure function 
at low region [9]. We denote this contribution as gY^^{x,Q^). 

(ii) In CSS collinear factorization scheme the soft gluons connect with the hard- and 
soft-parts can be absorbed into the soft-part at the collinear approximation, where the 
transverse momentum kr of the partons is neglected. However, the /cT-effects of the parton 
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at low should be considered. Thus, the collinear factorization should be replaced by 
the /cT-factorization scheme. Unfortunately, we haven’t a satisfy fcr-factorization scheme 
for the spin structure functions. The /cr-effects also include the replacements of Sq{x, Q^) 
and the DGLAP equation with the transverse momentum dependent (TMD) distribntion 
Sq{x, kTi Q^) and corresponding new evolution equations. While we haven’t such tools yet. 
We assume that a satisfactory choice of the parameters in the input parton distributions 
can mimic these /cj’-effects. 

(iii) According to the operator product expansion (OPE) in QCD [10], the Q^-evolntion 

of structnre function moments can be described in terms of twist expansion. The twist-2 
represents the scattering from individnal partons, while higher twist corrections appear 
dne to correlations among partons. At low scale, the higher twist (HT) contribntions 
to the structnre functions play a signihcant role. A special twist-4 corrections-the parton 
recombination to the DGLAP evolution eqnation at LL{Q‘^) approximation has been 
derived by us for un-polarized and polarized parton distributions in [4,5]. We will detail 
its contribution qDGLap+zrs ^ext section. While the typical contributions of the higher 
twist power suppressions ~ to gi are neglected in this work since 

they mainly change at a; > 0.1 [11]. 

(iv) The more complicated corrections to gi at low are from the higher order 
QGD effects 0(af). In our works we only consider the contribntions of the leading order 
corrections. An nnavoidable qnestion is whether we can neglect all higher order QGD 
corrections when <C IGeU^? Since all order resummation of these corrections are 
difficult, we take following generalized leading order (GLO) approximation [12]: if the 
leading order contribntions (or including necessary lowest order corrections) to a given 
process are compatible with the experimental data, one can conjecture that these neglected 
higher order corrections to this process may cancelable each other, or they are snccessfully 
absorbed by a hnite number of free parameters. 

In consequence, at small x and low we have 

g,(x, Q'-‘) ~ gf Q^) + QA (2,7) 
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We will detail every term of Eq. (2.7) in next sections. 
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3 Contributions of parton distributions 


In this section, we present the contributions of the polarized parton distributions of 
the proton to the spin structure functions at low i.e., 

gDaLAP+ZRS^^^g-2^ = lT.4Sg{x,Q^) + Q% ( 8 ) 

^ <? 

where the QCD evolution dynamics take the DGLAP equation with the parton recombi¬ 
nation (ZRS) corrections at the LL{Q‘^) approximation (see Eqs. (2.1)-(2.11)in Ref.[12]), 
the minimum free parameters in the input parton distributions have been fixed by the 
data mainly at a; > 10“^, they are 



xuy{x,jx^) = 24.3a;^'®®(l — 

(3.1) 


xdy{x,jx^) = 9.10a;^'^^(l — 

(3.2) 

in [13], and 




6uy{x, p^) = 40.3a;^'®^(l — a;)^'^^. 

(3.3) 


6dy{x,jj,‘^) = —18.22a;^'^^(l — x)*'^, 

(3.4) 

at p^ = 0.064GeR^ in 

[12]. The input distributions Eqs. (3.1-3.4) neglect the contribu- 


tions of asymmetry sea quark corrections. 

Note that 1/p = 0.78/m consists with a typical proton scale 0.8 — 1 fm. We consider 
that p is a minimum transverse momentum of the partons in the proton due to the 
uncertainty principle. Thus, we assume that all parton distributions are freezed at scale 
if < p^. Based on this assumption we avoid the un-physical singularities at 
Q ~ ^QCD- 

We present x-dependence of {x, Q‘^) at several values of in Fig. 2. One 

can End the dramatic change of the spin structure function at a; < 10“^ from a flat form 



to dramatically decreasing. Considering Fig. 3 in Ref. [12], we conclude that the large 
gluon helicity effect leads to this phenomenon. 
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4 Contributions of the VMD part 

As we have emphasized that the contribution from the vector meson in virtual photon 
to Qi at < IGeV"^ is necessary. Traditionally, this correction of the vector meson can 
be described by the VMD model and it had been used for the predictions of structure 
function at small x and low region. 

This contribution is written as 




mlQ^ 


- ^Cyp ( 5 ) , 


where 7 ^, is the coupling constant of vector meson and proton; a; is a variable dehned as 
X = Q‘^/{s + — rrip) rather than a momentum fraction of parton, s is the CMS energy 

square for the jp collision. The contributions of u meson are similar to that of p, while 
the contributions of 0 meson are small at < 1 GeV^ and can be neglected. We take 
V = p and a; at < IGeV^ and = 0.770GeV . The cross-sections Aa^pis) is the 
total cross section for the scattering of polarized meson with the nucleon, unfortunately, 
they are unknown. Usually, the following parameterized formula is used. 


A(Jyp(^S^ ^,at X < Xo 


(4.2) 


Thus, we take 


ln"^a; 


g=) . ,4.3) 

where the second factor aX x > xq is an arbitral function to suppress the contributions 
of the VMD mechanism with increasing x. The extrapolation of Pi from the measured 
region down to a; ~ 0 suggest us to assume that A = 1 — e and Xq = 10“^, where e ~ 0 is 
a small positive parameter due to the requirement of integrability of pf at a; —>• 0. In this 
work, we temporarily take e = 0. Thus, 
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VMD 


{x,Q^) = B 


mlQ‘^ 


+ mlY X 


[{^)0{xo -x) + ~ ^ 0 )], 


In Xo 


(4.4) 


where the parameter B = 3. 
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5 Predictions for spin structure function at small 

X 

What is the asymptotic behavior of This is a broadly discussed subject. We plot 
gi{x, Q"^) = Q‘^)+gY’^^{x, with different values of in Fig. 3. There 

are two different asymptotic behaviors of gi at x < 10“^: the VMD behavior ~ x~^ at 

< IGeV^^ and the large gluon helicity effect at > 3GeV^. Besides, gi presents the 
twist form of the two asymptomatic behaviors above, which is the mixing result of the 
nonperturbative and perturbative dynamics. 

We compare our predicted gf at x > 10“^ with the data [14] in Fig. 4. These data 
on 2010 are more precise than the previous data. Note that the values of of every 
measured point are different and they are taken from Table I of [14]. The theoretical curve 
is a smooth connection among these points. This figure shows that the pQCD evolution 
almost control the behavior of gi at x > 10“^. 

On the other hand, the combination of non-perturbative and perturbative dynamics 
at a; < 10“^ leads to a dramatic change of gl around = 1 ~ 3GeV^. Unfortunately, 
there are only several data with large uncertainty about gi in this range. In Figs. 5 and 

6 we collect the HERA early data [15,16] at = 1, lOGeU^ which are un-generally used 
and compare them with our predicted gf. Figure 7 shows some of these data (trigon) 
[15] and the comparisons with our results (dark points). Figure 8 is the Q^-dependence 
of gi with hxed x, the data are taken from [17]. One can End that our predicted gi are 
compatible with these data, although more precise measurements are necessary. 

Finally, we compare our results with the new COMAPSS (primary) data [1,2] at 

< lGeV‘^, which show that gi presents a flat asymptomatic form at a; < 10“^. This 
seems to contradict with the predicted strong rise of gl at < IGeV"^ in Fig.3. However, 
in the COMPASS fixed target experiments there is a strong correlation between x and 
Q^, which makes low x measurements also with low Q^. In Fig. 9 we take the average 
values of for each probing values of x (see Fig.l in Ref.[1]). The results are acceptable. 
Obviously, the measurements at different x with different values of in the fixed target 
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experiments mix two different asymptomatic behaviors of gi. 

We predict the stronger Q^- and x-dependence of gl at 0.01 < < 3GeV^ and a; < 0.1 

due to the mixtion of nonperturbative vector meson interactions and the QCD evolution 
of the parton distributions in Fig. 3. For testing this prediction, the measurements of gf 
with hxed x or at low are necessary. The planning Electron-Ion Collider (EIC), 
for example, eRHIC [18] and EIC@HIAF [19] can probe a broad low Q'^ < IGeR^-range, 
where we can check the predicted behavior of gi at hxed x or Q^. 
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6 Discussions 


In general consideration, both the logs oil/x and are equally important at small x 
and low and one should sum the double logarithmic (DL) terms (os ln^(l/a;))", which 
predict the singular behavior gl ~ x~^ (A > 0). To this end, some special attentions are 
proposed [20]. For example, the double logarithmic terms are taken into account via a 
suitable kernel of the evolution equations in the infrared evolution equations, which was 
hrst suggested by Lipatov [21], or alternatively taking a singular initial parton distribu¬ 
tions at a; < 10“^, one can also mimic the results of the DL-resummation. 

In this work, the behavior of Qi at the same range is obtained through a long evolu¬ 
tion of the DGLAP equation with the parton recombination corrections. We hnd that it 
is different from the predictions of the DL-resummation, the asymptomatic behavior of 
the polarized quark distributions at a; —)■ 0 is controlled by ^Pqg in the DGLAP equa¬ 
tion, rather than the ln^(l/a;)-corrections to the DGLAP-kernel. Thus, the difficult DL 
resummation can be replaced by the hts of the initial quark distributions Sqv{x,n‘^) in 
the DGLAP equation if the evolution distance is long enough. This conclusion was also 
obtained in the unpolarized structure functions [22]. 

In summary, we use the DGLAP equation with the parton recombination corrections 
and the nonperturbative vector meson dominance model to predict the spin structure 
functions Qi of the proton. We hrst present a complete picture for the translation of gf 
from low 0) to high at small x. We hnd that the contribution of the large gluon 

helicity dominates gi at x > 10“^, but the mixtion with nonperturbative component 
complicates the asymptomatic behavior of gi at x < 10“^. The results are compatible 
with the data including the early HERA estimations and GOMPASS new results. The 
predicted strong Q^- and x-dependence of gi at 0.01 < < ?>GeV^ and a; < 0.1 due to 

the mixtion of nonperturbative vector meson interactions and the QGD evolution of the 
parton distributions can be checked on the next Electron-Ion Gollider (EIG). 
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Figure 1: The time ordered decomposing of DIS diagrams, (a) The struck quarks are 
on-mass-shell since they have only forward component, (b) A ’’cat ear” diagram, which 
vanishes in the collinear factorization schema, (c) The ’’cat ear” diagram with higher 
order QCD corrections, which are non-vanished at low Q^, but can be isolated using a 
naive VMD model. 
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Figure 2: Perturbative contributions to g^. All partons are evolved from 

three valence quarks at = 0.064Ge\^^. 
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Figure 6: Predicted gf at = lOGeV^^ and the comparison with the HERA ’’data”, 
which are based on the NLO QCD predictions with the statistical errors expected at 
HERA [16], 
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Figure 7: Predicted g{ at 5.6GeV^^ and 16.5GeV^^ at a; < 10 ^ (circles) 

and the comparison with the HERA data (triangles) [15]. 
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Figure 8: Predicted Q^-dependence of with fixed values of x. Data are taken from [17] 





Figure 9: Predicted as a function of x with different measured Q‘^(x) (solid curve). 
Note that the low values of x connect with the low values of Q‘^{x). The data are taken 
from COMPASS primary results with two different beam energies [1]. 


24 




